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Probing Gaussianity represents one of the key questions in modern cosmology, because it allows 
to discriminate between different models of inflation. We test for large-scale non-Gaussianities 
in the cosmic microwave background (CMB) in a model-independent way. To this end, so- 

fj called first and second order surrogates are generated by first shuffling the Fourier phases 

belonging to the scales not of interest and then shuffling the remaining phases for the length 

scales under study. Using scaling indices as test statistics we find highly significant signatures 

for both non-Gaussianities and asymmetries on large scales for the WMAP data of the CMB. 

^H We find remarkably similar results when analyzing different ILC-maps based on the WMAP 

^ five and seven year data. Such features being independent from the map-making procedure 

would disfavor the fundamental principle of isotropy as well as canonical single-field slow-roll 
inflation - unless there is some undiscovered systematic error in the collection or reduction of 
the CMB data or yet unknown foreground contributions. 
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<^> 1 Introduction 

One of the key questions in cosmology is to probe the Gaussianity of the primordial density 
fluctuations, because it allows to discriminate between different models of inflation. While the 
simplest model of inflation, namely single- field slow-roll inflation 1 , predicts that the temperature 
fluctuations of the cosmic microwave background (CMB) correspond to a (nearly) Gaussian, 
homogeneous, and isotropic random field, more complex models may give rise to non-Gaussianity 
(NG). Models in which the Lagrangian is a general function of the inflaton and powers of its 
first derivative can lead to scale-dependent non-Gaussianities, if the sound speed varies during 
inflatioiEH Similarly, string theory models that give rise to large non-Gaussianity have a natural 
scale dependence^. Possible deviations from Gaussianity have been investigated in studies based 
on e.g. the WMAP data of the CMB (see e.g. Komatsu et alPand references therein) and claims 
for the detection of non-Gaussianities and other anomalies, like hemispherical power asymmetry, 
the 'axis of evil' , the Cold Spot etc. have been made. However, most of the tests on non- 
Gaussianities do not take into account any possible scale-dependency. In this contribution we 
apply the formalism of surrogate data sets, which was recently adapted to CMB data analysis 
9 , to test for non-Gaussianities on large scales in the WMAP data. Special emphasis is put on 
a comparison of the WMAP five year data with the WMAP seven year data. 
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Figure 1: Schematic view of the shuffling procedure on the /-ra-plane for generating the first and second order 

surrogates. The gray arrows indicate the phase permutations for obtaining first order surrogates and the black 

arrows show the shuffling for making the second order surrogates. 



2 Data Sets 

We use the foreground-cleaned Internal Linear Combination (ILC) maps generated and provided 
by the WMAP-team on the basis of the five year (WMAP5JP and seven year (WMAP7jr data. 
For comparison we also analysed the five year ILC-map (NILC5) produced by Delabrouille 
et alP, which was generated pursuing a needlet-based approach for removing the foreground 
contributions. 



3 Surrogates and Scaling Indices 

The model-independent test for scale-dependent non-Gaussianities is based on the use of so- 
called surrogate data sets. As test statistics for NGs we use scaling indices. Both methods were 
introduced and described in detail previously! 9 ! 10 ! 11 ! Here, we briefly review the main points of 
the two formalisms. 



3. 1 Surrogates 

Consider a CMB map T(0, 0), where T(0, <j)) is Gaussian distributed, which can easily be achieved 
by a rank-ordered remapping of the temperatures onto a Gaussian distribution. We calculate the 
Fourier transform of T(9,(j)). The complex valued coefficients a/ m , ai m = /dfi n T(n) 1^(71) can 



be written as a\ m — \ai m \e z ^ lrn with 



arctan(Jra(a/ m )/i2e(a/ m )). The linear or Gaussian 



properties of the underlying random field are contained in the absolute values |a/ m |, whereas all 
higher order correlations (HOCs) - if present - are encoded in the phases <^ m and the correlations 
among them. To ensure that the distribution of the phases is uniform, we also perform a rank 
ordered remapping of the phases (j)i m . To test for scale-dependent NGs we first generate a first 
order surrogate map, in which any phase correlations for the scales, which are not of interest 
(here: the small scales), are randomized. This is achieved by a random shuffle of the phases 
(f)l m for / > l cut , < m < /, where l cut = 20 throughout this study and by performing an inverse 
Fourier transformation. Second, N (here: N = 500) realizations of second order surrogate 
maps are generated for the first order surrogate map, in which the remaining phases (f)i m with 
1 < / < Z cu t?0 < m < I are shuffled while the already randomized phases for the small scales 
are preserved. In Fig. [I] the phase shuffling procedure for generating first and second order 
surrogates is schematically visualized. 





Figure 2: Deviation S as derived from rotated upper hemispheres for (a(rio)) for the WMAP5 ILC map (left), 
the WMAP7 ILC map (middle) and the needlet-based ILC map (right). 



3.2 Scaling Indices 

To compare the two classes of surrogates, we calculate local statistics in the spatial domain, 
namely scaling indices (SIM) as described previously 10 n . In brief, scaling indices estimate 
local scaling properties of a point set P. The spherical CMB data can be represented as a 
three-dimensional point distribution P — pi = (x^, y^, Zi), i = 1, . . . , N p i xe i s by transforming the 
temperature fluctuations into a radial jitter. For each point pi the local weighted cumulative 



point distribution p is calculated p(pl,r) = Yl 



3=1 



-m 2 



with dij = \\pl —pj\\. The weighted 



scaling indices a(pl,r) are then obtained by calculating the logarithmic derivative of p(pl,r) 
with respect to r, a(pi, r) — — ^wfr • Using the above-given expression for the local weighted 
cumulative point distribution p, the following analytical formula for the scaling index a 



a(pl,r) = 



j=i z y~r 
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is obtained. For each pixel we calculated scaling indices for ten different scales, ri 
0.25 in the notation of Rath et alP. 
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4 Results 



For each scale we calculate the mean ((a)) of the scaling indices a(pl,r) derived from a set of 
pixels belonging to rotated hemispheres. The differences of the two classes of surrogates are then 
quantified by the a-normalised deviation S(Y) = (Y surrol - (Y surro2 ))/a Ysurro2 , Y = (a(pl,rj)) : 
j = 1,...,10, surrol: first order surrogate, surro2: second order surrogate. Fig. [2] shows the 
deviation S((a(rio))) as derived from pixels belonging to the respective upper hemispheres for 
768 rotated reference frames for the three ILC-maps under study. The z-axis of the respective 
rotated reference frame pierces the center of the respective colour-coded pixel. 

Statistically significant signatures for non-Gaussianity and ecliptic hemispherical asymme- 
tries become immediately obvious, whereby the patterns of asymmetry remain remarkably 
similar for the three maps. Interestingly enough, we obtain slightly larger deviations for the 
WMAP7 map (S m in = -3.99, Sm ax = 3.73) as compared to the WMAP5 map (S m in = -3.87, 
Smax = 3.51). We find the largest deviation (S m i n = —4.36, S max = 4.5) for the NILC map, 
which can be considered as a more precise full-sky CMB temperature map than the ILC maps 
generated by the WMAP team 8 . Thus, the level of non-Gaussianity systematically increases, 
when the underlying CMB-map becomes less noisy and less foreground contaminated. 
To quantify the similarity of the patterns of asymmetry in the three maps we calculate the 
cross-correlation C(r) of the S(Y)-ma,ps derived from the three input maps as a function of the 
scaling range r. The results are shown in fig. [5] (left). For each scaling range r the three maps 
are highly correlated among each other with C(r) always being larger than 0.87. The highest 
correlations are found for the largest scales r, where C(r) reaches values of 0.98 and more. Thus, 
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Figure 3: Left side: Cross correlation coefficient C(r) for the 5 (y) -maps as a function of the scaling range r. 
Black: NILC vs. WMAP5, blue: NILC vs. WMAP7 and red: WMAP5 vs. WMAP7. Right side: Minimal and 
maximal values of the cr-normalized deviations S(Y) for the rotated hemispheres for each scaling range r. Black: 

NILC, blue: WMAP5 and red: WMAP7. 

on all scales r the patterns of asymmetry are very similar for each input map. On the right 
hand side of fig. [3] we show the minimum and maximum of S as a function of r. Except for the 
smallest r's we obtain for each map stable 3cr-deviations for both extrema, where - once again 
- the NILC map always yields the largest deviations for scaling ranges r > 0.1. 

5 Conclusions 

In conclusion, we detect highly significant signatures for asymmetries and non-Gaussianities 
for large scales (/ < 20) in the WMAP five and seven year data. The increase of the signal 
with decreasing noise as well as the very high correlations between the significance maps points 
towards an intrinisic nature of the detected anomalies, which are independent of the map mak- 
ing procedure. Such features would disfavor the fundamental principle of isotropy as well as 
canonical single-field slow-roll inflation - unless there is some undiscovered systematic error in 
the collection or reduction of the CMB data or yet unknown foreground contributions. Thus, 
further tests are required to rule out other systematic effects as origin of the detected anomalies. 
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